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1. INTRODUCTION 


Usually, in the works оп the electronic structure of the disordered alloys one 
concentrates on calculating the configuration average of the local density of states. 
Nevertheless, sites in a random alloys are not equivalent due to random fluctuations 
in the local atomic environment. Additionally, alloys generally have atomic short- 
range—order and for that reason the average value of the site's occupation number 
can be different from the occupation number of any particular site. It seems, that 
the static fluctuations of electronic charge density in metallic binary alloys has been 
overlooked in literature on alloys properties and only in work by Ling and Gelatt [1] 
these problems were studied intensively. On the other hand, the similar problems 
(from mathematical point of view) i.e. the magnetic moments fluctuations and its 
concentration dependence in ferromagnetic alloys were studied in detail by Hamada 
and Miwa [2]. Also many experiments have studied the effect of the local atomic 
environment on the local magnetic moments (see References in Ref. [2]). Among 
the experimental probes available to check or investigate charge fluctuations in 
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metallic alloys are, for example, core-level energy shift using X-ray photoemission 
spectroscopy, Móssbauer isomer-shift measurements, the optical measurements of 
the defect states in the semiconductors and others. 

In a work by Ling апа Gelatt [1] the static charge fluctuations were investigated 
by means of three different calculational techniques. These techniques were a weak- 
scattering perturbation expansion, a dilute-alloy perturbation expansion and a 
coherent potential approximation — perturbation theory. We can see, that the 
results obtained from these calculations are valid only in a narrow range of some 
parameters and do not take into account the short-range-order. It should be 
noted that actual alloys generally have atomic short-range-order. For that reason 
it would be interesting and at the same time very important to investigate the 
influence of the atomic short-range-order on the static charge fluctuations in the 
binary metallic alloys. More exact calculations of the charge fluctuations can be 
also very useful in estimation of the broadening of a core-level linewidth observed 
in X-ray photoemission (see also for experimental evidence of a fluctuating charge 
state in cupric oxide extracted from the X-ray photoelectron spectroscop'c study 
is). 

The paper is organized as follows. In the next section we give the general 
expression for static charge fluctuations. In Sec. III we present the method needed 
for calculation of the alloy electronic characteristics with the short-range-order 
included. The concrete evaluations of the different alloy's electronic characteristics 
and charge fluctuations are performed in Sec. IV. 


2. THEORY 


2.1. Charge fluctuations 7 


In this section we present several basic relations concerning the charge fluctuations 
in alloys [1]. For each microscopic arrangement of a type-A and type-B atoms (here 
we consider the substitutional binary alloys), all information about the electronic 
structure is contained in the Green's function 


G(E) = (E - Н)! (1) 


where H is a model alloy's Hamiltonian. This Green's function depends on a 
specific arrangement of a constituent atoms. Having in hand С one can calculate 
the site-projected density of states according to the formula 

М(Е)--І Im Tr (>< GE") (2) 


where |: > is the Wannier state at 1-th site and E+ = E +16, = > 0. In order to 
calculate the site occupation number one have to integrate over energy the site— 
projected density of states 


ni = f "© ar (E)N (E) dE (3) 
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where np (E) is the Fermi distribution function. Usually, as a main quantity charac- | 
terizing the charge distribution in an alloy, the average value of the site occupation 
number is calculated 


Ep 
пА(8) = A < МЕ) »;-A(p) dE (4) 
-00 

where < ... >;—4 denotes a restricted configuration average, Le. a configuration- 
averaging process is performed only over those arrangement with A atom at the 
origin. Неге we confined ourselves to the case of zero absolute temperature, so'the 
chemical potential is replaced by the Fermi energy Er. Of course, am represents 
the charge distribution rather in a crude manner. For that reason we have to 
calculate another quantity which can better reflect the true charge distribution. 
In this paper we are going to investigate the mean-square fluctuations in the site 
occupation numbers and according to the formulas given in [1] we can write 


Aną = ((ni- < ni Pica) іт 


Ep Ep 
+ f dE І ФЕ ((N;(E)- < N;(E) »;-A) (МЕ) < МДЕ) >а) а 
(5) 


An, == dE dE" < ( ImTr |i ><:|G(E) – ImTr |: >< i|G*74(E))- 
oo -00 


-( ImTr |i >< :|G(E') — ImTr |i >< i|G=A(E')) э-а (6) 


where G*=4 (E) denotes the alloy Green's function in a case for which the i-th site 
is occupied by A-atom. 
2.2. Cluster-Bethe-lattice method 


In order to calculate the influence of the short-range-order on the mean-square 
charge fluctuations (MSCF) according to the formula (6) we have to calculate a 
restricted average values of the product of two Green's function. 

In a study of disordered systems a widely accepted concept is that of an effective 
„medium. More specifically, a homogeneous medium is chosen in such a way that 
the average behaviour of an ensemble of disordered samples is reproduced in the 
behaviour of one homogeneous'sample. A widely used example of such approach is 
the coherent potential approximation (CPA) [4], or its generalizations, for exam- 
ple, the homomorphic CPA [5], the self-consistent self-energy approximation [6], 
the molecular CPA [7], the travelling cluster approximation (see [8] and references 
therein), and others. Here we use thé cluster-Bethe-lattice method (CBL) соп- 
nected with a single site CPA, which allows us to calculate in a relatively simple 
way the alloy Green's function with short-range-order included. The Bethe-lattice 
approximation used, for example, for the calculation of the partition function of a 
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magnetic systems [9] and for which an exact analytic treatment is possible for sim- 
ple tight-binding Hamiltonians. It consists in substituting for the infinite periodic 
lattice an infinite system of connected atoms with the same coordination number 
as the lattice of interest, but without closed rings of bounds. This lattice has the 
property that the one-particle Green's function at a given site can be expressed in 
terms of the Green's function at the preceding site in the lattice. The Bethe-lattice 
approximation is exact in one-dimensional systems and is good in low coordination 
three-dimensional lattices. For a closed-packed structures one have use a better 
approximation, usually the Husimi cactus model which consists of an infinite sy- 
stem of connected tetrahedra [10, 11]. In the CBL approximation a small cluster of 
atoms is treated exactly and the remainder of the alloy is simulated by attaching 
Bethe-lattice of the same coordination number *o the "dangling bounds" of the 
atoms at the surface of the cluster. The alloy Bethe-lattice can be, in a next step, 
replaced by the effective Bethe-lattice with the effective fields calculated within the 
single site approximation. 

Аз was already mentioned, we shall use the CBL method for calculation of the 
Green's function conditional average values appearing in Eq. (6). Let us consider 
a cluster of nine atoms (a smallest cluster of atoms for bcc lattice) which is descri- 
bed by the single-particle tight-binding (single-band) Hamiltonian with a hopping 
integrals V;; taking values V44, Увв, VAB, and Vga depending on whether the 
i-th and j-th atoms are both of class A, both of class B or one of each class, 
respectively. The short-range-order is introduced through the independent order 
parameter Рав (additional to the concentration of the one of a constituent atoms). 
This parameter denotes the possibility of finding a type-B atom as a'nearest ne- 
ighbour to a type-A atom. Six parameters, z, y = 1 — z, Рад; Рав, Ppa ала 
Рвв describing the concentration of the type-A and type-B atoms and the degree. 
of correlation in the occupation of a neighbouring sites can be expressed (because 
of symmetry and normalization conditions) by two independent parameters. № 
is most convenient to choose these to be z and p, where p is introduced by the 
following relations (see, for comparison, [2, 12 - 17]) 


Paa = 2+(1- 2)р 

Рав = (1—х)(1—р) (7) 
Psa = zc(1-p) 

Ppp = 1- z(1 - г) 4 


Note, that completely random disorder occurs when the short-range-order pa- 
rameter p is equal zero. For p — —1, z — 0.5 perfect long-range-order occurs and 
for p = 1 the alloy separates into pure- A and pure-B phases (segregation limit). It 
should be noted that the local order is extended only over a small region depending 
upon both z and p and usually extends to the fifth nearest-neighbours [18]. 
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In the Hamiltonian describing the Bethe-lattice alloy we have taken the average 
hopping matrix element in the form [14] 


У = z(PAA 'Vaa + Рав Van) + (1- z)(Pap Vas + Pea: Vea) (8) 


and the diagonal self-energy о(Е) was introduced to replace the random diagonal 
elements by site independent value [4]. Following the transfer matrix method (see, 
for comparison, [19]) and coherent potential approximation [4] the self-energy can 
be calculated from following equations: 


«(ел — g(E))(1 — (ел — о(Е))бог) "+ 


(1 — z)(es — o(E))(1 — (ев — о(Е)) Соо) = 0 (9) 
г Goo = 26% + (1 — z)G$o (10) 
СА = (E-e4 - ZYT)”! * (11) 
GB, = (E- eg - ZVT)”! (11) 


where the transfer matrix T' reads as 


p Ё 8) + У(Е - (Е)? -32V (13) 
14V 


Here, Goo is a diagonal matrix element of the average Green's function at site 0” 


and G$, is a diagonal matrix element of the average Green's function at site 20 


when this site is occupied by type-A atom. 

In the next step we consider a cluster of atoms (in our case Z + 1 atoms, Z = 8) 
and to each of the surface atom in this cluster a Bethe-lattice is connected. From 
the Dyson equation for a Green's function we have [14] 


Z 
(Е — =о) Со = 1 + Y VoiGio (14) 

:=1 
(E — ei)Gio = VioGoo + (2 – 1)Ӯ, КЧ, i=1,...,Z (15) 
(Е – с(Е)) Kl) =VGxv+(Z-1)VKf", i=1,...,2 (16) 
(Е – о(Е)) KP =УК + (2- ЦУК®), i=1,...,Z (17) 


where ки! is a Green's function connecting the :-th atom in the j-th level of the 
Bethe lattice with the central atom, V; = Уд ог Ув depending on a type of atom 
occupying the i-th lattice site and 


Va = РААИАА + РАВУАВ 
(18) 


Ув РввУвв + РвАУвл 


Il 
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Having in hand equations (7 - 18) and the expression for the average value of a 
Green's function 


2 
8 n —Пл пА.2 -ЋА . 
ТО АК (w, ) С енеден qu) 9 


і-0 


. we are able to calculate the average values needed for calculation of the charge 
fluctuations. Eq. (6). Here, GAZA Б) denotes the diagonal matrix element 
of a Greens function G, Eq. (1), in a case when in the central site of a cluster 
is located a type-A atom surrounded by n4 type-A atoms and (Z — n4) type- 
B atoms (first neighbours). Additionally, with each of these boundary atoms the 
effective Bethe-lattice is connected. The explicit expression for Gata) (Е) is 
as follows (cf. [14, 15]): 


бал"; Е; ke - ea) A "АУДА - (8— па)Уль. —. (20) 
2 , 
X=(E- ex) РА, (21) 
E — o(E) — Ж 
Үр фу. М (22) 
| E - e(E) - — : 
_ E-o(B) + V (E - o(E))? - 280” (23) 


2V 
To the end, the expression for the mean-square fluctuations in the site occupation 
numbers reads as [1]:. 


А 1 f Ep Ер ; 
An (в) = = dE | dE' < Im Goo(E) : Im Goo(E") >o= в) — 
T \ ж T 


Er 1 
-(f dE < Im Ggo(£) 2272 , (24) 


~ where the calculations of the conditional averages are performed according to the 
formula (19). 


3. NUMERICAL RESULTS AND DISCUSSION 


In this section we apply the method of Sec. 2 to study the charge fluctuation in 
an alloy with the diagonal, off-diagonal disorder and short-range-order included. 
We present also the density of states corresponding to a single type-A and type- BH 
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atoms located inside of a nine atoms cluster and surrounded by the effective Bethe- 
lattice. In order to better visualize the set of a short-range-order parameters p and 
a type-A atoms concentrations z, for which the calculations have been performed, 
we present its as a small circles in Fig. 1 in the region of (z, p) parameters (see for 
comparation, [12]). As one can see, we choose such values of the short-range-order 


0 0.125 0255 05 P> 1 


иди 


Fig. 1. Allowed range of the short-range-order parameter p vs. concentration of the type-A 
atoms х. The small circles denote points for which calculations of the charge fluctuations were 
done 


parameters in order to investigate interesting for our study phenomena at various 
limiting cases, ie. near the segregation limit (for constant value of z), at random 
disorder (for increasing values of z), for constant value of z and increasing values 
of p and in additional positions near the perfect order limit. In our calculations we 
have taken Ед = —ep = 2.0, 1.0 or 0.5 eV and for Vaa, Vas = Vsa and for Увв 
we have chosen the set of values (1.0, 1.1, 1.2) eV (see, for comparison [13, 18]). 
In all cases presented here we have only one energy band (without localized states 
outside of the band limits). 

In Figs. (2-6) we have drawn the local density of states (LDOS) corresponding 
to the central type-A and type-B atoms of the nine-atoms cluster-Bethe-lattice 
system for all possible distributions of the atoms in a cluster. Additionally, we 
have presented also the average LDOS — the broken line. From Figs. (2, 3) one 
can observe the changes of LDOS for equal concentrations of the constituent atoms 
and for disorder parameter going from the value corresponding to the long-range- 
order to random disorder case. Note that in Fig. 2. one has P44 = Рвв = 0, 
Рав = Pga = 1 and in Fig. 3, Pag = Рав = Pga = Рвв = 0.5. In Fig. 2, for 
type-B atom in a centrum of the cluster constructed from type-B atoms only, we 
find a pronounced peak in lower part of the energy band. The strength of this peak 
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Fig. 2. Local density of states corresponding to the central type- A atom (left-hand side) and 
type- В atom (right-hand side) in the nine-atoms cluster surrounded by effective Bethe-lattice. 


The number of different atoms in the cluster is indicated as (na, ng) for na = Ù, 1,...,8, 
where ng = 2— ПА, Z = 8. The broken curve denotes the average local density of states. 
The concentration of type-A atoms т = 0.5 and short-range-order parameter р = —1 


gradually decreases as we move from long-range-order limit towards the random 
disorder. Also, the average LDOS changes significantly. For a long-range-order we 
have symmetrical LDOS which changes into rather structureless curve for random 
disordered alloys. Looking on Figs. (3 — 5), one can observe a continuous changes 
of LDOS with increasing value of a type-A atom concentration. 

Here we have Pag = Рав = Pga = Pap = 0.5, z = 0.5 for Fig. (3), Paa = 
Pga = 0.75, Рвв = Рав = 0.25, z = 0.75 for Fig. (4) and P44 = Pga = 0.825, 
Pap = Рвв = 0.125, z = 0.875 for Fig. (5). There is a noticeable change (with 
increasing value of z) of the LDOS curves corresponding to all types of the clusters. 
The average LDOS has a two-peaked structure in the vicinity of the binary-ordered 
limit. 

In Fig. (6) we have shown the same characteristics but for z = 0.25 near 
the segregation limit, ie. Paa = 0.8125, Рьв = 0.9375, Рав = 0.1875 and 
Рвл = 0.0625. In general, we observe most peaked structures in LDOS in a 
case of random disorder and more smooth LDOS curves when some value of the 
short-range-order is introduced. Having calculated LDOS for different clusters we 
proceed to computing the mean-square charge fluctuations in the site occupation 
numbers. Results are shown in Figs. (7 - 11). 


We display the charge fluctuations versus the Fermi energy level. In Figs. (7 — 8) 
we have displayed MSCF for z — 0.25 and increasing value of a disorder parameter 
on type~A atoms (upper panel), here minority sites, and on type- B atoms (lower 
panel), here majority sites. The site energies are c4 = -ер = 0.5 eV and 2.0 eV for 
figures (7) and (8), respectively. We observe the largest fluctuations for a case with 
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Fig. 3. The same as in Fig. (2) but for p = O (random disorder limit 


largest disorder, ie. for points which lie relatively close to the random disorder 
limit. We also observe, that for a greater scattering strength (SS) the MSCF are 
greater than for smaller values of SS. In addition, for small value of SS the MSCF 
has a form of the two peaked structure, without third maximum in upper part of 
the energy band. For greater SS the largest fluctuations appear in the middle and 
upper part of the band. In Fig. 9 we present MSCF for points on (x,p) plane 
(see Fig. (1)) lying along the line of a constant concentration of a type-A atoms, 
т = 0.5 and varying disorder parameter р (from nearly long-range-order limit to 
a random disorder case). For a point z = 0.5, p = —1 we have obtained MSCF а 
few order smaller than this one for greater values of the parameter p and we do not 
displayed it in tnis figure. The greatest MMSCF are observed for р = 0 and when 
we move towards smaller values of the parameter p (p « 0) then the fluctuavions 
are decreased. Only in a lower part of the energy band (for the fluctuations on 
a type-B atoms) one can observe another sequence of peaks (in comparison with 
MSCF on a type-A atom) but differences are rather small In Figs. (10, 11) we 


1 А 5 ; 
present MSCF for p — rE and т = 0.5, 1 = 0.75. Тһе site energy c4 is 1.0 eV in 


Fig. (10) and 2.0 eV in Fig. 11. First of all, we observe the general feature that 
with increasing value of the scattering strength the charge fluctuations increase. 
Additionally, for smaller values of SS and far away from the random disorder limit, 
the MSCF has a structure with two maxima (and not with three, as usual), without 
the maximum placed for Fermi energies lying in higher part of the energy band. 
In conclusion, we have presented a theory of the static charge fluctuations in 
alloys with a short-range-order included. .Тһе charge fluctuations, as well as the 
local density of states was calculated within the method which treats part of the 
system exactly as à cluster of atoms and simulates the rest of a lattice by con- 
necting to each cluster's surface atom a Bethe-lattice. Generally, the mean-square 
charge fluctuations calculated for a system with the short-range-order included are 


134 у В. Taranko, E. Taranko, M. Pilat 


-6-4-2 0 2 4 6 8-6 -4 -2 0 2 4 6 
E lev] Elev) 


Fig. 4. The same as in Fig. (2) but for z = 0.75 and р = 0 


reduced in comparison with fluctuations obtained for completely random disorder 
' alloy. The results presented in this paper indicate, that the charge fluctuations for 
a long-range-order limit are a few orders smaller than for other values of the para- 
meter describing the alloy's disorder. Also, for points (z, p) lying near the binary 
ordered limit the fluctuations are enhanced on majority sites, whereas on minority 
sites are depressed. Finally, the influence of the short-range-order on the charge 
fluctuations in alloys is significant, so any theoretical description of the experimen- 
tal study of the real alloys systems for which the amount of a valence charge on an 
atom can be important, have to take this fact into account. 
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Fig. 5. The same as in Fig. (2) but for z = 0.875 and p = 0 


Fig. 6. The same as in Fig. (2) but for z = 0.25 and p = 0.75 


136 В. Taranko, Е. Taranko, M. Pilat 


Fig. 7. Meen-square charge fluctuations on type-A atoms (upper pane?) and on іуре-В atoms 


(lower panel) as a function of Fermi energy calculated for three different values of the 


short-range-order parameter р = 0.5, 0.75 and 0.875 (curves A, В and C, respectively) for 
2 = 0,25 and £4 = —ep = 0.5 eV 


7:8 © -2 0 2 & Ew) 


Fig. 8. The same as in Fig. (7) but for £4 = —ep = 2.0 eV 
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; D 1 
Fig. 9. The same as in Fig. (7) but Юг z = 0.5 and p = 3 3 and 0.0 (curves A, B and 
С, respective") for £4 = —ЕВ = 1.0 еу 
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| 


Fig.10. The same as in Fig. (7) but for p — -% 1 = 0.5 апа 0.75 (сигуез А, апа В. 


respectively) for £4 = —Eg = 0.5 eV. The charge fluctuations on minority atoms (i.e. оп 
type- DB atoms, lower panel) are a few orders smaller and not displayed here 


. 
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Fig. 11. The same as in Fig. (10) but for £4 = —€g = 1.0 eV 
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STRESZCZENIE 
t 

W pracy przedstawiono teorię fluktuacji ładunkowych w stopach dwuskładnikowych, chara- 
kteryzujących się istnieniem uporządkowania bliskiego zasięgu. Uporządkowanie to uwzględniono 
przy pomocy dokładnego rozważenia małej grupy atomów, natomiast wpływ pozostałej części kry- 
ształu opisany jest efektywną siecią krystaliczną typu drzewa Cayley'a. Przeprowadzono obliczenia 
numeryczne lokalnej gęstości stanów dla klasterów dziewięcioatomowych (o różnej konfiguracji 
atomów) oraz fluktuacji ładunku dla szerokiego zakresu parametrów charakteryzujących stopy. 
Fluktuacje ładunku są największe w przypadku stopów, które nie wykazują uporządkowania bli- 
skiego zasięgu i znacznie maleją w miarę jego wzrostu. Duże fluktuacje ładunku zaobserwowano 
również w stopach wykazujących duże uporządkowanie typu А-В-А-В 

W świetle uzyskanych wyników należy sądzic, że teoretyczny opis zjawisk fizycznych, w 
których wielkość ładunku zlokalizowanego na atomach stopu ma zasadnicze znaczenie, musi 


uwzględniać wpływ uporządkowania bliskiego zasięgu na wielkość tego ładunku. 


РЕЗЮМЕ 


В работе представлено теорию флуктуации заряда в бинарных сплавах харак- 
теризующихся существованием ближнего порядка. Упорядочение это учтено с 
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помощью точного описания малой группы атомов, а влияние остальной части 
кристалла описано эффективной кристаллической решёткой типа дерева Кейлея. 
Проведено машинные расчеты локальной плотности состояний для 9-ти aTOM- 
ных кластеров (о различной конфигурации атомов) и флуктуации заряда в случае 
широкого предела параметров описывающих сплавы. Флуктуации заряда самые 
большие для сплавов, в которых не существует ближний порядок и значительно 
уменьшаются по мере его возрастания. Большие флуктуации наблюдают 
PM цыдот ся тоже в сплавах, которые характеризуются большим упорядочением 
типа А-В-А-В. 

Из наших результатов следует, что теоретическое описание физических яв- 
лений, в которых величина эффективного заряда локализованного на атомах 
сплава играет существенную роль, должно учитивать влияние ближнего поры- 


адка на величину этого заряда. 


Zlozono 13. I. 1989. 
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